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ABSTRACT 

In the article a test is developed, which allows to test the null-hypothesis of 
the intrinsic randomness in the angular distribution of gamma-ray bursts col- 
lected at the Current BATSE Catalog. The method is a modified version of 
the well-known counts-in-cells test, and fully eliminates the non-uniform sky- 
exposure function of BATSE instrument. Applying this method to the case of 
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all gamma-ray bursts no intrinsic non-randomness was found. The test also 
did not find intrinsic non-randomnesses for the short and long gamma-ray 
bursts, respectively On the other hand, using the method to the new in- 
termediate subclass of gamma-ray bursts, the null-hypothesis of the intrinsic 
randomness for 181 intermediate gamma-ray bursts is rejected on the 96.4% 
confidence level. Taking 92 dimmer bursts from this subclass itself, we ob- 
tain the surprising result: This "dim" subclass of the intermediate subclass 
has an intrinsic non-randomness on the 99.3% confidence level. On the other 
hand, the 89 "bright" GRBs show no intrinsic non-randomness. 
Subject headings: cosmology: observations - gamma rays: bursts 



1 INTRODUCTION 

Two results of the last years in the statistics of the gamma-ray bursts (GRBs) 
are doubtlessly remarkable. 

The first one concerns the number of subclasses. Recently, two differ- 
ent articles ( ||Mukherjee et al. 1998|| ; |[Horvath 1998|1 ) simultaneously suggest 



that the earlier separation |[Kouveliotou et al. 1993|| of GRBs into short and 
long subclasses is incomplete. (It is a common practice to call GRBs having 
T 90 < 2 s (T 90 > 2 s) as short (long) GRBs, where T 90 is the time during 



which 90% of the fluence is accumulated |[Kouveliotou et al. 1993|| .) These 



articles show that, in essence, the earlier long subclass alone should be fur- 
ther separated into a new "intermediate" subclass (2 s < T90 < 10 s) and 
into a "truncated long" subclass (T 90 > 10 s). (In what follows, the long 
subclass will contain only the GRBs with T90 > 10 s, and the intermediate 
subclass will be considered as a new subclass.) 

The second result concerns the angular distribution of GRBs. At the 
last years several attempts ( [[Hartmann et al. 1991|| , fBriggs et al. 199"GH ; 



legmark et al. 1996b|| ; [Balazs et al. 1998]; [[Balazs et al. 1999|j ) were done 



either to confirm or to reject the randomness in the angular sky distri 
bution of GRBs being collected at BATSE Catalog ( ||Fishman et al. 1994| |: 



Mccgan et al. 1998|| ). Theoretically, if the intrinsic distribution of GRBs is 



actually random, an observation of some non-randomness is still expected due 
to the BATSE non-uniform sky-exposure function ( ||Fishman et al. 1994 



Meegan et al. 1998H ). Hartmann et al. (1991), Briggs et al. (1996) and 



Tegmark et al. (1996b) did not find any statistically significant depar- 
ture from the randomness. On the other hand, the existence of some non- 
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randomness was confirmed on the > 99.9% confidence level by Balazs et al. 
(1998). This behavior can be caused either purely by instrumental effects 
or the instrumental effects alone do not explain fully the detected behavior 
and some intrinsic non-randomnesses should also exist. Balazs et al. (1998, 
1999) suggest the second possibility. This conclusion follows from the result 
that while the short subclass shows a non-randomness, the intermediate + 
long subclasses do not indicate it. It is difficult to explain such behavior of 
subclasses by the instrumental effects alone. 

In this article we will again investigate the angular distribution of GRBs. 
Trivially, after the discovery of the new intermediate subclass, it is highly 
required to test the intrinsic randomness in the angular distribution of this 
new subclass, too. In addition, of course, new different tests, which exactly 
eliminate the effect of the sky-exposure function, are also required in order 
to complete the results of Balazs et al. (1998, 1999). 

The aim of this article is to test the intrinsic randomness in the angu- 
lar distribution of all GRBs and of the three subclasses separately, too. We 
will use a modification of the well-known counts-in-cells method. This is a 
standard and simple statistical test (see, cf., Meszaros (1997) and references 
therein). The advantage of this method is given by the fact that it allows 
to eliminate quite simply and exactly the sky-exposure function. The main 
result of paper will be the surprising conclusion that the intermediate sub- 
class and only this subclass alone suggests a non-randomness on the 96.4% 
confidence level; its "dimmer" half even on the 99.3% confidence level. 

The paper is organized as follows. In Section 2 the method is described. In 
Section 3 the results of test are presented. Section 4 discusses and summarizes 
the results of the article. 

2 THE TEST 

Assume for the moment that there is no non-uniform sky-exposure function. 
We separate the sky in declination into m dec > 1 stripes having the same area 
(47T /m dec steradian). The boundaries of stripes are the declinations 5 k, k — 
0, 1, mdec, where 8 = —90 degree and 5 mdec = +90 degree, respectively. 
The remaining values are analytically calculable, and appear symmetrically 
with respect to 5 = 0. One has: sin 5k = 2k/md ec — 1- (For example: If 
m dec — 3, then 5 12 = ±19.47 degree; if = 4, then 81^ = ±30.00 degree 
and 5 2 = 0.00 degree; etc.) We also separate the sky in right ascension a 
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into m ra > 1 stripes. They are defined by boundaries a = 360/c' /m ra degree; 
k' = 0, 1, ..,m ra . (Obviously, a trivial modification of this separation in right 
ascension is the case, when the boundaries are a = 360 +p)/m ra degree, 
where p is an arbitrary real number fulfilling < p < 1 .) All this means that 
we separated the sky into M = rridec x mra areas ("cells") having the same 
size Ait /M steradian. If there are N GRBs on the sky, then n = N/M is the 
mean of GRBs at a cell. Let rif, i = 1,2, ...,M be the observed number of 
GRBs at the zth cell (E^i m = N). Then 



Al 



var M = (M-l)- 1 J2(ni-n) 2 (1) 



defines the observed variance. For the given cell structure with M cells, 
due to the Bernoulli distribution ( ||Meszaros 1997|| , [Balazs et al. 1998fl ), the 



measured variance varu should be identical to the theoretically expected 
value n(l — 1/M). This theoretical prediction should then be tested. 

Note that this and similar methods (see, e.g., Meszaros (1997) for details 
and further references) are usual in astronomy. For example, this method was 
used already by Abell (1958) to reject the randomness in the sky-distribution 
of clusters of galaxies. The test is, compared with other statistical tests 
("two-point angular correlation function", "nearest neighbor distances", etc.; 
[[Peebles 1980|| ; fDiggle 1983|| ; [Pasztor 1993|| ), not the most sensitive one to 



detect non-randomnesses. Its importance for our purposes is given by the fact 
that it allows an extremely simple generalization to the case with non-zero 
sky-exposure function. 

Now, we generalize the method to this case. This may easily be done 
by changing the boundaries of cells in order to have the same probability 
(and hence the same expected number n = N/M) for a given cell. The sky- 
exposure function is a function of declination only ( ||Fishman et al. 1994j| ; 



Meegan et al. 1998|| ). Hence, the choice of equatorial coordinates is highly 



convenient, because then no changes of boundaries are necessary in right 
ascension. The new boundaries 5k, k = 0, 1, ...,md ec in declination may be 
calculated analytically as follows. Clearly, 5o,m dec = ±90 degree remain. In 
BATSE Catalog [|Meegan et al. 1998| the exposure function f(5) is defined 



for 37 values of declination (for 5 r = — 90, — 85, +85, +90 degree; r = 
0, 1, 2, .., 36). To obtain 5k we, first, calculate the value 

A^—J2f(5 r )cos5 r , (2) 

i8U r=l 
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where for the given r the corresponding declination is 5 r = —90 + 5r degree. 
(Remark that r = and r = 36, respectively, need not be in the sum, 
because cos(±90) =0.) Then, second, for rridec > 2 we search for the values 
5k] k = 1, 2, (rridec — 1) as follows. For the given k we search for the 
declination Si fulfilling the condition 

7 E/(^)cos^<-^<-J 7 EM)cos5 J . (3) 



36 A £1 rn dec 36 A ^ 

Having this we search, by linear interpolation between Si and for the 
exact value of 5k- By this method 5k is well calculable. (For example: For 
m dec = 3 we obtain 5i = —19.51 degree, 5 2 = 22.44 degree; for rridec = 4 we 
obtain S\ = —30.83 degree, 5 2 = 1.51 degree, 5 3 = 33.60 degree; etc.) Having 
these cells with these " shifted" boundaries in declination the variance may be 
calculated identically to the case with no sky-exposure function. This method 
will test the pure intrinsic randomness; the effect of BATSE sky-exposure 
function is exactly eliminated. 

It is natural to probe different values of M. In addition, for some M dif- 
ferent cell structures are still possible (cf. M — 12 allows m dec = 2,3,4,6). 
Hence, generally, several - say Q - cell structures may be probed for the same 
sample of GRBs. Having these Q cell structures (and hence Q means + Q 
measured variances) two questions arise. 1. How to calculate the confidence 
level for a given cell structure? 2. Having Q values of confidence levels, how to 
calculate the final confidence level? The answer for the first question seems 
to be quite clear: varuln seems to be identical to the y 2 value for M — 
1 degree of freedom ( [|Trumpler fc Weaver 19531 ; jKendall fc Stuart 1969 



Press et al. 1992| ; the mean is obtained from the sample itself, and there- 
fore the degree of freedom is M - 1). Nevertheless, the situation is not 
so obvious, because the x 2 test needs n > 5 ( flrumpler fc Weaver 1953 



Kendall fc Stuart 1969|| ; ||Press et al. 1992|| ). In addition, some statistical 



text-books propose to use "quadratic" cells only Qpiggle 1983|| , Chapt. 2.5.). 
If all these restrictions were taken into account, then \ 2 tests would be pos- 
sible only for 2m rfec = m ra ; M = 2m dec and iV > 5M = 10m dec . This would 
be a drastic truncation of the possible cell structures. But, not doing these 
restrictions, the estimation of the confidence level for a given cell structure 
must be done by more complicated procedures; e.g., by numerical simula- 
tions. Concerning the answer to the second question the situation is even 
less clear. As the reasonable search for the final confidence level only Monte 



Carlo simulations seem to be usable IPress et al. 1992 
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Keeping all this in mind, we will proceed as follows. In the coordinate 
system with axes x — 1/M versus y = J votm / n = (var/mean) 1 ^ 2 the Q 
values of (var jmean) 1 ! 2 define Q points (one point for any cell structure; 
Vj = \JvarM,j/n, where j = 1, 2, Q). Clearly, for these points one expects 
the theoretical curve y = a/1 — x. This theoretical expectation can straight- 
forwardly be verified, e.g., by least squares estimation ( [Press et al. 1992| , 



Chapt. 15.2.; [piggle 1983|| , Chapt. 5.3.1). Our estimator is the dispersion 



a Q = £(%' - V 1 " l / M ) 2 • (4) 
i=i 

Obviously, smaller ctq suggests that the theoretical curve is better fitted. 
Note still that, as the best choice, the square root of var M /n is proposed in 
this "var/mean" test ( ||Diggle 1983] , Chapt. 5.4.). 

The confidence level can then be estimated by Monte Carlo simulations 
in the following way: We throw 1000-times randomly N points on the sphere, 
and repeat the above calculation leading to oq for every simulated sample. 
Then we compare the size of the oq obtained from this simulation with <jq 
obtained from the actual GRB positions. Let u be the number of simulations, 
when the obtained <tq is bigger than the actual value of oq. Then one may 
conclude that (100 — u;/10) is the confidence level in percentage. Clearly, this 
method does not need n > 5 and quadratic cells. 

There is no commonly accepted confidence level in statistics, above which 
the null-hypothesis should already be rejected ( ||Trumpler fc Weaver 1953 



Kendall fc Stuart 1969|| ). It is only a general agreement that confidence lev- 



els smaller than 95% should not be considered. Our opinion is (see also 
Kendall fc Stuart 1969|| ) that the confidence levels bigger than 95% can al- 



ready be taken as "remarkable", "suspicious", "interesting", etc.; a higher 
than 99% confidence level may still mean the rejection of null-hypothesis, 
and such result must doubtlessly be announced. Hence, we will require that 
the confidence level be bigger than 95%. Thus, here it must be u < 50. 

In this paper GRBs will be taken between trigger-numbers 0105 and 6963 
from Current BATSE Catalog [[Meegan et al. 1998|| having defined T 90 (i.e. 



all GRBs detected up to August 1996 having measured T 90 ). From them 
we exclude, similarly to Pendleton et al. (1997) and Balazs et al. (1998), 
the faintest GRBs having a peak flux (on 256 ms trigger) smaller than 0.65 
photon/(cm 2 s). This truncation is proposed by Pendleton et al. (1997) in 
order to avoid the problems with the changing threshold. The 1284 GRBs 
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obtained in this way define the "all" class. From them there were 339 GRBs 
with T 90 < 2 s (the "short" subclass), 181 GRBs with 2 s < T 90 < 10 s 
(the "intermediate" subclass) and 764 GRBs with T 90 > 10 s (the "long" 
subclass). We will study the all class and the three subclasses separately. 

We will ad hoc choose rridec = 2, 3, . . . , 8 and m ra = 2, 3, . . . , 16. I.e. it will 
be Q = 105. Of course, this choice of Q is more or less subjective. Neverthe- 
less, our choice is motivated by two concrete arguments. First, we would like 
to study only the angular scales much bigger than the positional errors. (The 
size of a cell will not be smaller than 22.5 degree. On these angular scales no 
problems should arise from the positional errors ||Meegan et al. 1998|| .) Sec- 
ond, it is reasonable not to consider such high values of m^ec, when 180/md ec 
is already comparable or even smaller than 5 degree. (If this were not re- 
quired then the elimination of sky-exposure function would be problematic 
due to its definition for declination intervals with widths 5 degree.) 



3 THE RESULTS 

Figure 1 collects the results of Q = 105 "var/mean" tests of four different 
cases. It is obvious immediately that for the "all" case the points follow 
well the theoretical curve. For the "short" and "long" subclasses, on the 
other hand, there is a slight tendency of points to be above the theoreti- 
cal curve. The situation concerning the intermediate subclass seems to be 
the most unambiguous: mainly for small M (roughly below M ~ 40) the 
points are clearly above the theoretical curve. This suggests an intrinsic 
non-randomness mainly in the sky distribution of intermediate subclass; such 
possibility for the short and long subclasses, respectively, cannot be excluded, 
too. 

The results of Monte Carlo simulations support this expectation only in 
the case of intermediate subclass. We obtain u = 287 (u = 80, uj = 36, uj = 
440) for all GRBs (short, intermediate, long GRBs). Hence, the rejection of 
null-hypothesis is confirmed for the intermediate subclass only on the 96.4% 
confidence level. For the short and long subclasses, respectively, and also 
for all GRBs the null-hypothesis cannot be rejected on the > 95% confidence 
level. For the short subclass we have a 92% confidence level; for the remaining 
two cases even smaller levels. 
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4 DISCUSSION AND CONCLUSION 

The most surprising result of paper concerns the intermediate subclass. The 
intrinsic non-randomness is confirmed on the confidence level > 95%. This 
confidence level, as discussed in Section 2, is "remarkable", but is not enough 
to reject the null-hypothesis of randomness. 

The results concerning the 339 short GRBs should also be mentioned. 
Nevertheless, the 92% confidence level is clearly not enough to reject the 
confidence null-hypothesis. On the other hand, this result, together with 
Balazs et al. (1998, 1999), suggest that also for the short subclass itself the 
rejection of null-hypothesis of intrinsic randomness can also occur by further 
tests. 

In the case of 764 long GRBs, and also of the 1284 all GRBs, there are 
no indications for the non-randomnesses. All this seems to be in accordance 
with the results of Balazs et al. (1998, 1999). 

We think that the result concerning the intermediate subclass is highly 
surprising, because just this new subclass, having the smallest number of 
GRBs, has a remarkable "proper" behavior. 

A short further investigation of this subclass fully supports this conclu- 
sion. 

There are 181 GRBs in this intermediate subclass. Be divided this sub- 
class into two further subclasses; into the "dim" and "bright" ones. By 
chance the peak flux = 2 photons/(cm 2 s) (on 0.256s trigger) is practically 
identical to the medium of peak flux for this subclass. Therefore, we consider 
the GRBs having smaller (bigger) peak flux 2 photons/ (cm 2 s) as the "dim" 
("bright") subclass of the intermediate subclass. There are 92 GRBs at the 
"dim" subclass, and 89 GRBs at the "bright" one. 

We provide the 105 "var/mean" tests for these two parts, too. We obtain 
the surprising result that the "dim" subclass has an intrinsic non-randomness 
on the 99.3% confidence level (u; = 7). Contrary this, the "bright" subclass 
can still be random (u = 662). 

The sky distribution of 92 intermediate dim GRBs is shown on Figure 2. 

We mean that the behavior of the intermerdiate subclass of GRBs, quite 
independently, supports the correctness of the introduction of this new sub- 
class ( |^ukher]ee^F^irT99^| ; RHorvath 199S|| ). Further investigations of this 



new subclass are highly required. 
Three notes are still needed. 

First, purely from the statistical point of view, it must be precised that 
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even the rejection of null-hypothesis of the intrinsic randomness would not 
mean a pure intrinsic non-randomness in the spatial angular distribution of 
GRBs. This is given by the fact that, up to now, it cannot be fully excluded 
that GRBs (or some part of them) are not unique phenomenons, and there 
can occur some repetitions, too. This question is studied intensively by 
several papers ( ||Mcegan et al. 1995 ], | Quasmiock 1995 ], | Quashnock 1996 ], 



Tcgmark et al. 1996a|| , ||Graziani et al. 1998j| , ||Hakkila et al. 1998|| ) conclu- 



ding that repetition can still play a role. 

Second, strictly speaking, the statistical counts-in-cells test is testing the 
"complete spatial randomness" (shortly the "randomness") of the distribu- 
tion of GRB on the celestial sphere ( IPiggle 1983 1 , Chapt.1.3). Therefore, in 
this paper we have kept this terminology. In cosmology, on the other hand, 
the word "random" ("non-random") is rarely used, and the word "isotropic" 
("anisotropic") is usual (for the exact definition of isotropy in cosmology see, 
e.g., Weinberg (1972), Chapt. 14.1). Of course, here we will not go into the 
details of these terminology questions (see, e.g., Peebles (1980) for more de- 
tails concerning these questions). We note only that the " random- isotropic" 
("non-random-anisotropic") substitution is quite acceptable on the biggest 
angular scales; on smaller angular scales the situation is not so clear. There- 
fore, in Balazs et al. (1998, 1999), where only the angular scales ~ 90 degrees 
and higher were studied, the words "isotropy" and "anisotropy" were quite 
usable. In this article, going down up to the scales ~ (20 — 25) degree, the 
used terminology is more relevant. 

Third, trivially, further studies are needed. They should test - by other 
different statistical methods - again the intrinsic randomnesses (more gener- 
ally: the intrinsic spatial distributions |[Lamb 1997]| ), both for all GRBs and 
for the subclasses. In addition, a test of the repetition alone, i.e. a test not 
being influenced by positions, is highly required. 

As the conclusion, the results of paper may be summarized as follows. 



• We developed a method, which can verify quite simply the intrinsic 
randomness alone in the angular distribution of GRBs, because the 
method eliminates exactly the non-zero sky-exposure function. 

• We rejected the null-hypothesis of the intrinsic randomness in the an- 
gular distribution of 181 intermediate GRBs on the 96.4% confidence 
level. 
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• We rejected the null-hypothesis of the intrinsic randomness in the angu- 
lar distribution of 92 "dim" intermediate GRBs on the 99.3% confidence 
level. 

• We did not reject the null-hypotheses of the intrinsic randomnesses in 
the angular distribution of the remaining two subclasses and of the 
all GRBs, respectively, on the > 95% confidence levels; the "bright" 
intermediate GRBs seem to be distributed randomly, too. 
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Figure 1. The results of 105 "var/mean" tests of four different cases drawn 
in the 1/M vs. \J var / mean frame. The theoretical curve sfl - l/M (solid 
line) is also shown. M is the number of cells. 
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Figure 2. Sky distribution of 92 GRBs of "dim" subclass of the intermediate 
subclass in equatorial coordinates. 
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